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Simulations of the Effects of Stripping and Accretion on 
Galaxy Haloes in Clusters. 
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ABSTRACT 

We present results from a series of hydro dynamic simulations investigating ram pres- 
sure stripping of galactic haloes as the host galaxy falls radially into a cluster. We 
perform a parameter study comprising of variations in initial gas content, gas injection 
rate (via stellar mass loss processes), galaxy mass and amplitude of infall. From the 
simulation results we track variations in both physical quantities (e.g. gas mass) and 
directly observable quantities (e.g. X-ray luminosities). The luminosity of the galaxy's 
X-ray halo is found to compare favourably with the observationally determined cor- 
relation with optical blue band luminosity (Lx : Lb) relation. Factors affecting the 
X-ray luminosity are explored and it is found that the gas injection rate is a domi- 
nant factor in determining the integrated luminosity. Observational properties of the 
material stripped from the galaxy, which forms an X-ray wake, are investigated and 
it is found that wakes are most visible around galaxies with a substantial initial gas 
content, during their first passage though the cluster. We define a statistical skewness 
measure which may be used to determine the direction of motion of a galaxy using X- 
ray observations. Structures formed in these simulations are similar to the cold fronts 
seen in observation of cluster mergers where a sharp increase in surface brightness is 
accompanied by a transition to a cooler region. 

Key words: galaxies: clusters: kinematics and dynamics - galaxies: interactions - 
intergalactic medium - X-rays: galaxies - galaxies: clusters: Virgo 



1 INTRODUCTION 

Most galaxie s are not isol ated but reside in clusters or groups 
of galaxies (iTullvl Il987l) where a number of mechanisms 
may affect the nature and evolution of the galaxy. These 
mechanisms include tidal interactions, galaxy harassment, 
and ram pressure stripping. Ram pressure stripping, which 
is the subject of this paper, has been identified in X-ray 
observations of several nearby elliptical galaxies e.g. M86 



jRangaraian et alJ Il995l IWhite et alJ 


19911 iForman et alJ 


I197SD and NGC 4472 Jlrwin & Sarazin 


Il99d) in the Virgo 



2002). In these relatively nearby cases it is possible to di- 
rectly image the X-ray emission surrounding the galaxy, 
which is distorted into an X-ray wake by the stripping pro- 
cess. Spectral analysis can yield information on the temper- 
ature and metallicity structure lending further weight to the 
hypothesis that the wake is the re sult of a stripping proce ss 
rather than an accretion process jRangaraian et allll99^ . 

To detect the effects of ram pressure stripping in more 
distant galaxies requires a different approach. As wakes are 
often low surface brightness features it is difficult to image 



wakes in distant systems, particularly when viewed against 
a bright cluster background. If stripping affects global prop- 
erties, such as the total X-ray luminosity of the galaxy, 
this could be detected, for example in an environmental de- 
pendence of the Lb relation. It is currently unclear how a 
galaxy's X-ray luminosity Lx is affected by the cluster envi- 
ronment and, if so, whether ram pressure is a significant or 
dominant mechanism. Different studies of elliptical galaxies 
in clusters do not always draw the same conc lusions (e.g. 
lO'Sullivan et alJl200ltlBrown fc Bregmanl feoOOi. It is possi- 
ble that the nearby wake candidates are atypical and that 
ram pressure stripping is not a common process. Alterna- 
tively, it may be the case that Lx suppression is not a good 
indicator of the presence of stripping or that other processes 
affect Lx more than ram pressure. 

Numerical simulations of stripping processes provide 
valuable insights into physical processes and their observ- 
able consequences which cannot be obtained from analytical 
models alone. Previous work in this area includes simula- 
tions of elliptical galaxies moving through a uniform cluster 
environment at a constant velocity which is ap propriate for 
galaxies on a circular orbit within a cluster fe.g. lGaetz et alJ 
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ll987l : lBalsara et al.lll994T) . Given the constant ram pressure 
in these cases a quasi-steady state is reached and certain 
scaling r elations (e.g. for re tained gas mass fraction) may be 
derived jGaetz et alJll987l) . Although Ikent fc Gurml (Il982l) 
determine that galaxy orbits in the Coma cluster cannot 
be primarily radial, the case of a radial, or highly ellipti- 
cal orbit, is still of considerable interest in representing ac- 
cretion of material from the field or a sub-cluster merger. 
NGC4839 is the dominant galaxy in a sub-cluster which 
is fallin g into the Coma cluster and XMM-Newton obser- 
vations ijNeumm^^^d] [gOO^, iNeumann et al] l2001) show 
an X-ray tail around NGC4839 which is interpreted as a 
wake of stripped material. As clusters grow via minor merg- 
ers there may be many such cases of large el liptical g alax- 
ies pa ssing through clusters on radial paths. Tak eda et alJ 
Jl984h simulate radial infall for a large (1O 12 M ) ellipti- 
cal galaxy tracking the galactic gas evolution for number 
of cluster crossings, finding a periodic variation in gas con- 
tent after the first cluster passage. A detailed study of el- 
lip tical galaxies on elliptical orbits has been carried out 
bv lToniazzo fc Schindleri j200ll) using 3-dimensional simula- 
tions yielding detailed information on the gas dynamic pro- 
cesses and variations in X-ray luminosity. A broad parameter 
stu dy of stripping from dw arf galaxies has been carried out 
bv iMori fc Burkerd ( 2000) for galaxy masses in the range 
10 6 -10 10 Mq . Interactions between dwarf galaxies and the 
cluster environment affect the chemical evolution of galax- 
ies and clusters; higher metallicity gas which is stripped from 
the galaxy will enrich the lower metallicity medium of the 
cluster, and the depletion of the gas content may prevent 
further star formation. Ram pressure stri pping may als o in- 
fluence the morphology density relation (Dressier 1980). 

Our aim in the present study is to use 2-D hydrody- 
namic simulations of a spherical galaxy falling radially into 
a cluster in order to determine the effect of such an event on 
the galaxy and the impact on observable properties such as 
the galaxy's X-ray luminosity. A number of parameter vari- 
ations are investigated enabling determination of the sensi- 
tivity of the results to the value of a given model parame- 
ter. This parameter study indicates not only how sensitive 
our results are to the choice of model, but also what are 
the consequences of genuine variations in the physical pa- 
rameters of a galaxy population. The aim of the parameter 
study presented here is to focus on high mass galaxies where 
interactions with the intra-cluster medium lead to more di- 
rectly observ able consequences than in the lower mass cases 
considered bv lMori fc Burkertl i2000t) . X-ray observations of 
galaxy wakes caused by stripped material could be used to 
elucid ate the full 3-D nature of galaxy dynamics iMerrifieldl 
1998) but a sample of wakes observed in a given cluster will 
not be unbiased hence it is important to understand the 
factors which influence which galaxies will host observable 
wakes. 

Section |21 presents the details of the galaxy model, the 
model cluster into which the galaxy falls and the numeri- 
cal scheme used. Section |3] consists of a number of subsec- 
tions which present the results from the simulation param- 
eter study. Subsection 13.11 gives an overview of the strip- 
ping phases which occur during a cluster crossing, subsec- 
tion examines the variations of the galactic gas content, 
the observable consequences of which are dealt with in sub- 
section 13.31 Subsection 13.41 examines the properties of the 



wake of stripped material and subsection 13.51 explores the 
properties of the bow shock. The effect of placing the model 
galaxy in the harsher stripping environment of a hotter clus- 
ter is examined in section^] Finally, section|K]presents a dis- 
cussion of the results and conclusions which may be drawn. 



2 OVERVIEW OF THE MODEL 

This section presents details of the canonical, or typical, 
simulation run. The parameter study involves varying each 
parameter from its value in the canonical model in order to 
determine the effect of altering a single model component. 
For example, to investigate the effect of increasing the size 
of the pre-existing gas halo, the initial extent of the halo was 
doubled whilst all other parameters retained their standard 
values. When changing the mass of the model galaxy, more 
than one simulation parameter was altered. The mass of 
the dark matter halo was changed and other parameters 
calculated so as to be consistent with the increased dark 
matter mass as described below. 



2.1 The galaxy model 

The model galaxy consists of two distinct components; a 
stellar component and a dark matter component. The dark 
matter component is considered to be the more massive and 
hence dominates the gravitational influence of the galaxy, 
however the stellar component is significant not only due 
to its gravitational effects but also due to sources of gas 
injection which come from stellar mass loss processes (e.g. 
type la supernovae, planetary nebulae and stellar winds). 
These injection processes place additional gas onto the grid 
within the galaxy and can replenish a stripped gas halo. A 
model galaxy is constructed by specifying a mass for the 
dark matter halo and then using a number of scaling rela- 
tions to determine values for subsequent parameters. 



2.1.1 Dark matter component 

Numerical s i mulati ons of structure formation by 
iNavarro et alJ Jl997t) indicate that dark matter haloes 
have a universal cuspy density profile where the dark 
matter density at a distance r from the centre of the halo 
is given by 

Po 



Pdn 



where 



cx (1 + cxY 



^200 



(1) 



(2) 



c is the concentration parameter which governs the degree 
of central concentration of the halo. In hierarchical models 
of structure formation smaller objects form earlier in the 
evolution of the universe and hence are more centrally con- 
centrated so the concentration parameter scales with halo 
mass. The scaling relation used was t hat determined from 
the simulations of iBullock et alJ \2QQ ll) . 
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The value of T2oo depends on the redshift of formation of the 
object and the matter density of the Universe at that red- 
shift which in turn depends on the assumed cosmology. Due 
to the difficulties in determining these quantities a typical 
value of T2oo = 50 kpc was chosen for Mhaio = 2.0 x 10 12 Mq 

1/3 

and it was assumed that r2oo oc M halo . The final constraint 
imposed on the dark matter halo is that the dark matter 
distribution is truncated at r2oo, so that 

1-^200 

M halo = / p dm (r) 47rr 2 dr. (4) 
Jo 

2.1.2 Stellar component 

Given the properties of the dark matter halo it is possible to 
define a stellar population which is consistent with the ob- 
served properties of elliptical galaxies. The total blue band 
luminosity of the stellar population was determined assum- 
ing a total mass to light ratio of 5 Mq /Lq based on Table 1 
from iFinoguenov fc Jonesl fcOOttl (N.B. this is the mass to 
light ratio for the whole galaxy, values measured using an 
aperture smaller than r2oo will measure smaller mass to light 
ratios as dark matter dominates at large radii). For the old 
stellar population of an elliptical galaxy the mass to light 
ratio of the stars themselves will be greater than unity. The 
value adopted was 



T7 



= 13- 



(5) 



based on an ave r age value from the sample of 
IFinoguenov fc Jonesl (l2000t) w here the orig inal mass to 
light ratios were determined bv lLauerl lll985l) . Given these 
assumed mass to light ratios the total stellar mass may 
then be determined. Two mass to light ratios are required 
as the first dictates the amount of dark matter relative to 
the galaxy's blue band luminosity and the second describes 
the mass to light ratio of the stars only which depends on 
the age of the stellar population. The stars are assumed 
to be d istributed according to a King-type distribution as 
used bv lBalsara et al.1 (119941) i.e. 



P* (r) 



P* (0) 



(6) 



1 + (r/RcY 

where R c is the core radius (2.5 kpc) and the stellar distribu- 
tion is truncated at Rh = 20 kpc. The central stellar density 
p* (0) is determined using 



Af* = 



p* (r) 47iT dr. 



(7) 



continuity between the outer edge of the halo and the sur- 
rounding intra-cluster medium (ICM). 



2.1.4 Mass replenishment and drop out 

Stellar mass loss processes (e.gtype la supernova (SNIa), 
planetary nebulae) will inject matter and energy into the 
galaxy which affects the X-ray halo. The distribution of re- 
plenished material is assumed to follow the stellar distribu- 
tion so that the mass injection rate as a function of radius 



Prep 



P* (0) 



1 + (r/Rc 



(8) 



where a* = 5.4 x 10 20 s 1 is the specific mass loss 
rate for an old stellar p opulation as determined by 
iBrighenti fc Mathews! J1998I) . 

The energy injection is dominated by SNIa and a 
component due to the velocity dispersion of the stars 
in the galaxy potential. A su pernova rate typical of an 
elliptical galaxy is 0. 1 SNU (IFinoguenov fc Jonesl l2000l : 
ICappellaro et al.1 Il993h (1 SNU = 1 supernova per year 
per Lb©) and the typical energy injection from one SNIa 
is 8 x 10 50 ergs. 

The velocity dispersion of the galaxy was estimated us- 
ing the Faber- Jackson relation 



2 x 10 10 L 



( Z V 

V200 km s" 1 / 



(9) 

200 km s _1 / 

flFaber fc Jacksonlll97rf l where Lv was determined from the 
blue band luminosit y Lb assuming a (B — V) colour of 0.95 
jBender et alJll993h . 

In addition to processes which act as sources of gas there 
are other processes which act as sinks to remove gas (e.g. star 
formation or accretion onto a central compact object). Ra- 
diative cooling is included, using a Raymond- Smith thermal 
plasma model with solar metallicity, which results in cool 
dense gas accumulating at the centre of the galaxy where 
it is allowed to drop out of the flow and is removed over 
a given time-scale. In a cell where the gas density is above 
the critical value of pdo = 2.4 x 10 _24 gcm -3 and the gas 
temperature is below the critical value of Td = 2 x 10 4 K 
gas will be removed at a rate of 

P 



Pdo = 



Tdo 



(10) 



where t^ = 3 x 10 7 y r. The drop out presc ription used here 
is identical to that of IStevens et alJ Ijlfl 



2.1.3 Gas halo 

The pre-existing gas halo, which surrounds the galaxy at 
the start of the simulation, is assumed to extend out 
to r2oo and be compos ed of material of solar metallicity 
llMatsushita et al]l2000l) . The halo is initially isothermal at 
a temperature of 0.8 keV, the measured temperatu re of the 
gas halo surrounding M86 jMatsushita et alJ I^OOO) which is 
of a similar mass to our model galaxy. The effects of radiative 
cooling are included which causes the dense central regions 
to cool appreciably from the initial temperature. The initial 
conditions are such that the halo is in hydrostatic equilib- 
rium with the galaxy potential and that there is pressure 



2.2 Cluster model and galaxy dynamics 

The model cluster was chosen to be a favourable environ- 
ment for observing wakes. In order for a significant number 
of wakes to be observed it is required that the density of the 
ICM be high enough that ram pressure is sufficient to distort 
the galaxy halo, but not so high that the background emis- 
sion from the cluster obliterates the signature of the wake. 
The optimum rang e of cluster tempera tures to fulfil these 
criteria is 2-3 keV dStevens et al]ll999ft . The parameters of 
the cluster model were based on observations of the cluster 
Abell 160, a dynamically relaxed cluster at a temperature 
of 2.7 keV and a distance of 270 Mpc, in which a statistical 
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Time (Gyr) Time (Gyr) 



(a) Time variation of galaxy position relative to the cluster 
centre. Positive values indicate that the galaxy is located to 
the right of the cluster centre and negative values indicate that 
the galaxy is located to the left. 



(b) Time variation of galaxy velocity. Horizontal lines mark the 
cluster sound speed of 847 km s — 1 



Figure 1. Galaxy dynamics. 



detec tion of a number of wakes has been found dPrake et alJ 
l200d) . 

The dynamics of the galaxy were modelled by free-fall 
in the cluster potential, neglecting dissipation processes such 
as dynamical friction but including the effect of tidal forces. 
The form of the cluster potential was chosen to be the same 
as that used by iToniazzo fc Schindlerl l|2C 



> (r) = ct c In 



1 + 



(0 



(ii) 



where <r c is the cluster velocity dispersion and r c is the clus- 
ter core radius. The free fall frame of reference in which the 
simulation is conducted transforms away first order gravi- 
tational effects but not tidal effects. For a potential of this 
form, tidal forces are greatest at the cluster core radius. A 
core radius of 400 kpc was adopted to match the value of 
IWhite et all l)l997h for Abell 160. The galaxy was allowed 
to fall into the cluster, on a radial trajectory, from an initial 
position of 2 cluster core radii from the centre and an ini- 
tial state of rest, resulting in the dynamics shown in Fig. Q 
The shape of the cluster potential is fixed by eqn. llll but the 
depth of the potential, parameterized by the cluster velocity 
dispersion, is not fixed. The X-ray surface brightness S(w) 
at a projected distance w from the centre of a relaxed cluster 
is found to be well fitted by a beta model with /3 = 2/3 



S(w) 



So 



(12) 



[l + (w/r c f 

which results from a gas density profile of the form 
Po 

1 + (r/r c ) 2 

Gas in hydrostatic equilibrium in the potential given in 
eqn, llll will have the density distribution given in eqn. 1131 if 



(13) 



kT 



(14) 



where fiiriu is the mean mass per ICM particle. A value of 
a c — 656 km s _1 was used, derived from eqn. 1141 with a tem- 
perature of 2.7 keV. The central ICM number density was 
determined by relating the total X-ray luminosity Lx and 
the central surface brightness So for the surface brightness 
profile given in eon. 1121 



Lx = 8n 2 Sor 2 . 



(15) 



The central surface brightness may be related to the central 
emissivity and hence to central number density. Assuming 
a plasma composed of hydrogen and helium in the ratio 
10:1 the total central number density was esti mated to be 
6.17 x 10" 4 cm" 3 for L x = 2.7 x 10 43 erg s _1 dWhite et alJ 
1997). 



2.3 Numerical scheme 

The numerical scheme chosen for this application was the 
VH-1 implementation of th e Piecewise-Parabolic Method 
JColella fc Woodward! Il984h which offers good resolution 
of shocks in complex flows and has b een extensively used 
in a variety of astrophys ical c ontexts jBk2ndjn_et L ^Jj|T990| , 
iRichards fc Ratlifj 1 19981 and IStrickland fc Stevens! boOOT) . 
The numerical grid was a uniform Cartesian grid of 800 x 400 
cells with a size of x = 1.6 x 10 24 cm, y — 8.0 x 10 24 cm where 
the centre of the galaxy was located at x = 8.0 x 10 24 , y — 0. 
In an axi-symmetric case, such as this, it is only necessary 
to simulate the top half of the galaxy and the symmetry 
of the situation can be used to reflect the flow field in the 
i-axis to obtain the full flow field which results in a signif- 
icant saving in computation time. The size of the grid as a 
whole was chosen such that features of interest, such as tails 
of stripped material and shocks, should not propagate off 
the grid. As the galaxy moves both from right to left and 
from left to right during the course of the simulation it is 
necessary to centre the galaxy with respect to the a;-axis. 
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The cells subtend an angle of ~ 0.5 arcsec at a distance of 
270 Mpc (the distance of Abell 160) which matches the spa- 
tial resolution of the Chandra X-ray observatory, currently 
the best available X-ray spatial resolution. 

Boundary conditions at the left and right of the grid 
were fixed inflow or fixed outflow conditions to account for 
the movement of the simulation frame in the gravitational 
potential of the model cluster (e.g. fixed inflow at the left 
boundary and fixed outflow at the right boundary for a 
galaxy moving from right to left). The magnitude of the 
inflow velocity at one boundary was equal to the magnitude 
of the outflow velocity at the other boundary so that the 
direction of motion could be reversed in a consistent fash- 
ion. The fixed inflow/outflow densities were different due to 
density gradients in the ICM. The size of the grid is not neg- 
ligible compared to the size of the cluster so the ICM den- 
sity is different at the left and right boundaries. To ensure a 
physically consistent model for the cluster, it is necessary to 
ensure the ICM contains a density gradient consistent with 
hydrostatic equlibrium in cluster potential given by eqn. llll 
If ean. H4l is obeyed then the resulting X-ray surface bright- 
ness distribution of the ICM will be given by a beta model 
with f3 — 2/3 as in eqn. 1121 This density gradient must be 
present in the initial ICM distribution and must also be ac- 
counted for when calculating inflow and outflow densities at 
the left and right boundaries. 

A reflective boundary condition must be applied to the 
y — boundary in order for the symmetry of the problem 
to be correctly applied and an inflow/outflow condition was 
applied to the top of the grid to allow material to flow freely 
off the grid should it reach the upper boundary. Time steps 
were determined by application of the Courant-Friedrichs- 
Lewy (CFL) condition with a CFL number of 0.8. 

The use of a 2D numerical grid results in a signifi- 
cantly shorter computation time than for a 3D grid, for 
comparable resolution, which is an important issue when 
performing the many runs required to cover a wide range 
of parameter space. The rate of stripping may depend 
on the formation of Kelvin-Helmholtz shear instabilities 
which can behave differently in 2D and 3D simulations. 
iBriiggen fc Hillebrandtl feOOlal) compare the mixing rate 
due to such instabilites in 2D and 3D simulations and find 
that the steady state mixing rate is very similar, support- 
ing the validity of the 2D treatment used here. It should 
be noted, however, that if magnetic field effects are included 
there is significant disagreement betw een 2D and 3D simula- 
tions teriiggen fc HillebranddboOlrJ) . As no magnetic fields 
are explicitly included in these simulations we proceed with 
a 2D model. 



3 RESULTS 

In this section results are presented from the parameter 
study examining both observational quantities, such as X- 
ray luminosity and surface brightness profiles, and other 
quantities (e.g. gas mass) which, whilst not directly observ- 
able, play an important role in determining the nature of the 
galaxy. The parameters involved in the study were: galaxy 
mass, mass injection rate, initial gas halo extent and the 
starting position of the galaxy (hence the maximum veloc- 
ity reached) as shown in Tabled 



3.1 Initial cluster crossing 

Density distributions from the canonical run are shown in 
Fig. |21 to illustrate various phases of the stripping process. 
Fig. |2(a)| is the initial gas distribution (i.e. a spherically sym- 
metric halo in hydrostatic support) which is beginning to be 
stripped in Fig. |2(b)| where the motion is subsonic. The halo 
presents a sharp edge in the upstream direction but strip- 
ping, via Kelvin-Helmholtz instabilities, prevents a sharp 
delineation in other regions. Fig. |2(c)| is taken from a su- 
personic regime and a bow shock, and associated density 
enhancement, can be seen. The final frame (Fig. |2(d)} is 
taken at the end of the first cluster crossing and density en- 
hancements are seen in front of the galaxy. As the galaxy 
approaches the edge of the cluster, and the ram pressure 
decreases, material which was previously confined behind 
the galaxy moves forward causing the enhancements seen in 
Fig.pfdH 

3.2 Gas content 

The variation in gas mass within the central 40 kpc, for 
galaxies with different initial halo extents, is plotted in 
Fig. |3(a)| The vertical lines mark regions where the galaxy 
velocity is supersonic (indicated by Sup.) close to core pas- 
sage. A more extensive gas halo is more massive not only due 
to its larger volume but also due to the higher central density 
which is necessary to satisfy the conditions of hydrostatic 
support and pressure continuity at the halo truncation ra- 
dius. Both the canonical and extended haloes experience an 
initial period during which the pre-existing halo is stripped 
before moving into a cyclic phase where the enclosed gas 
mass is at a maximum at core passage. The duration of 
the initial stripping phase depends on the halo mass at the 
start of the simulation (~ 2.5 Gyr for the canonical halo 
and ~ 4.0 Gyr for the larger halo) but after this phase all 
runs behave in a very similar fashion indicating that after 
~ 4 Gyr (time of second core passage) the signatures of the 
initial gas content are no longer present. The large enclosed 
mass close to core passage is somewhat surprising, but one 
should bear in mind that this gas is not necessarily grav- 
itationally bound to the galaxy. Fig. |4(a)| shows the mass 
of gas which is gravitationally bound to the galaxy for the 
three cases considered above. The gas in a cell where the 
mass density is p, the particle number density is n, the fluid 
velocity u and the temperature T is considered to be bound 
if 

ipu 2 + ^nkT + pifjo < (16) 

where ipc is the gravitational potential of the galaxy evalu- 
ated at the cell location. In contrast to Fig. |3(a)| which plots 
all gas within 40 kpc irrespective of whether it is bound, 
Fig. |4(a)| indicates that the supersonic phase results in a 
substantial amount of the bound gas being stripped. This 
indicates that around the time of core passage there is a 
concentration of gas within 40 kpc which is not gravitation- 
ally bound to the galaxy and the question of where this gas 
is located is raised. Fig. plots the ratio of gas density af- 
ter 3.9 and 3.6 Gyr, i.e. the gas density at one of the peaks 
in Fig. |3(a)| divided by the gas density just after the peak 
starts to form. This indicates that the rise in unbound gas 
mass within 40 kpc is associated with an enhancement of 
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Table 1. Details of runs forming the parameter study 



Parameter 


Low Value 


Standard Value 


High value 


Halo extent 


25 kpc 


50kpc 


100 kpc 


Mass replenishment rate 


0.12 M yr" 1 


1.2 Mq yr" 1 


3.6 M yr" 1 


Dark matter halo mass 


1.0 x 10 12 M Q 


2.0 X 10 12 M 


4.0 X 10 12 Mq 


Infall radius 


400kpc 


800 kpc 


1600 kpc 




X (kpc) X (kpc) 

(c) t = 1.6Gyr, d = -350kpc (d) ( = 2.4Gyr, d = -790kpc 

Figure 2. Gas density distribution from the initial infall of the canonical run. The time (t) quoted is that since the start of the simulation 
and d is the distance from the cluster centre which is positive when the galaxy is to the right of the cluster centre and negative when to 
the left. The logarithmic greyscale shows the gas density in units of gem -3 . 



material behind the shock front. When the halo extent is 
reduced by half, the initial period of stripping prior to the 
cyclic behaviour is not seen as the initial gas content does 
not exceed that present during the latter phases as in the 
other cases. 

The gas mass for runs with different mass replenish- 
ment rates is shown in Fig. |3(b)| In the previous series of 
runs the initial gas content was different but the injection 



rate was constant resulting in the evolution of the gas con- 
tent in the model galaxies converging. In this second case the 
initial gas content is constant but the injection rate varies, 
resulting in a gas content which displays an initial period of 
divergence. At later times (after ~ 3 Gyr) the gas content 
again behaves in a cyclic manner but with a base level which 
depends on replenishment rate. The rate at which the mass 
increases prior to the supersonic phase is a function of the 



Galaxy Haloes in Clusters 7 



" 2.0, 

O 

2 

l 

a 

2 1-8 h 



O 1.4 
ft 

5 ^ 



'E 



Sup. 



Sup. 




mp 



ft 



2 





Sup. 




Sup 


1 


Sup 




\ 


M 
// \ i 
y i \\ 

' V 














(: 

V 
\\ 
\N 
\ 








<v 

1 b 






\ 


-A 

\ \ / 

\ ^ y 

\ 

./ 


v ^ 

1 


i 


h \ 


^\ 

X _ .— 



(a) Halo extent 



Time (Gyr) 
(b) Replenishment rate 



" 1 .8 
o 

1.6 
1.4 

O 
ft 

M 1.2 

o 

1.0 
3 .6 

1 .6 
CO 

m -1 

c0 -4 

6 , 

cn 
cO 

O o 





1 

Sup. 


i i 


i 

Sup 


i 


Sup. 


i 


















/ \ 

\ 




//\\ 
/ \\ 




//\\ 
/ \\ 






\ 

i 


i i 


i 


i 


\ 





Time (Gyr) 
(c) Galaxy mass 




1 2 3 4 5 

Time (Gyr) 



6 7 8 



(d) Infall amplitude 



Figure 3. Time variation of the gas mass within 40 kpc of the galactic centre. In each case the canonical run is plotted as a solid line, 
the dashed line represents a reduction in the parameter and the dotted line an increase in the value of the parameter. Time periods 
where the galaxy velocity is supersonic are marked "Sup.". Periods of supersonic velocity are not shown in Fig. |3(d)] as they are different 
for each case. 
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Figure 5. Gas density field after 3.9 Gyr divided by gas density 
field after 3.6 Gyr showing the location of the density enhance- 
ment around the time of core passage. The galaxy is moving from 
left to right at both times. 



mass replenishment rate with higher injection rates result- 
ing in more rapid rises in gas content. The higher injection 
rate run shows how in subsonic phases there is an accumu- 
lation of gas, during the supersonic phase the gas content 
is high due to the accumulation of unbound gas, and after 
the supersonic phase much of the gas associated with the 
galaxy has been depleted and is replenished during the next 
subsonic phase. The bound gas content (Fig. |4(b)| better 
shows the effect of alternate cycles of net gas accumulation 
during subsonic phases and net stripping during supersonic 
and some subsonic regimes. When the mass injection rate 
is at its lowest it is seen that the net stripping phase starts 
sooner than in other runs. During subsonic motion, strip- 
ping does occur from the outer edge of the halo but the 
net gas budget is a combination of the stripping and drop 
out loss mechanisms, and the mass injection replenishment 
mechanism. Whether there is net loss or gain of gas during 
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Figure 4. Total gas mass which is gravitationally bound to the galaxy (c.f. Fig. EJwhich plots all gas within a fixed 40 kpc radius). In 
each case the canonical run is plotted as a solid line, the dashed line represents a reduction in the parameter and the dotted line an 
increase in the value of the parameter. 



a subsonic phase will depend largely on the effectiveness of 
subsonic stripping compared to the mass injection rate. 

The next series of runs involves variation of the galaxy 
mass which is achieved by changing the dark matter halo 
mass and assuming length scales vary as r oc M 1 ^ 3 . By vary- 
ing the total galaxy mass, not only is the depth of the po- 
tential well increased but also the mass replenishment rate 
(which is proportional to the stellar mass), the mass replen- 
ishment temperature (which depends on the stellar velocity 
dispersion) and the extent of the initial halo (which is sub- 
ject to the same linear scaling factor as the other galaxy 
length scales) . The resulting interdependency of parameters 
means that more care is required in interpreting the results 
of these particular simulation runs but gives a self consistent 
galaxy model. The gas mass fraction, within a 40 kpc radius, 
for each galaxy is plotted in Fig. [(J where the gas mass frac- 
tion has been calculated by dividing the gas mass within 
40 kpc by the total galaxy mass (dark matter and stellar 
components). The gas mass fraction indicates that the more 
massive galaxies have a higher gas content in the subsonic 



regime. The enhancement due to the bow shock, during su- 
personic phases, is greater for less massive galaxies. As the 
concentration parameter scales inversely with galaxy mass, 
the less massive galaxies will be more centrally concentrated 
and this may result in a more prominent shock. The bound 
gas content (Fig. |4(c)| indicates that more massive galax- 
ies not only have more gravitationally bound gas but also 
the variation in this gas content is less. The deeper gravita- 
tional potential wells of the more massive systems would be 
expected to make depletion of gas by ram pressure stripping 
less effective. 

The final parameter to be varied is the initial starting 
position from which the galaxy falls into the cluster. This af- 
fects the galaxy dynamics (i.e. the maximum speed reached 
by the galaxy and also the time taken to cross the cluster, 
hence the periodicity in mass variation) . When the infall ra- 
dius is reduced to 400 kpc the maximum velocity reached is 
770 kms -1 and the crossing time is ~ 1.7 Gyr; in this case 
the galaxy velocity never becomes supersonic. An increased 
infall radius of 1600 kpc raises the maximum velocity to 
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Figure 6. Gas mass fraction within 40 kpc of the galactic centre 
for galaxies with differing masses (solid: canonical, dashed: less 
massive, dotted: more massive). This is essentially Fig.[3]with the 
i/-axis divided by the total mass of the galaxy. 
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Figure 7. Time variation of the gas mass within 40 kpc of the 
galactic centre for galaxies falling into the cluster from different 
starting positions with the time axis scaled to the cluster crossing 
time (solid: canonical, dashed: infall from closer to the cluster 
centre, dotted: infall from further out from cluster centre). 



1560 kms -1 with a crossing time of ~ 4.5 Gyr. The initial 
state of the gas halo also varies in these runs. When the 
starting position of the galaxy is closer to the cluster core, 
where the ICM density is higher, the halo density must also 
be higher to maintain pressure continuity at the halo bound- 
ary and hydrostatic support. The time variation of the gas 
mass for these runs, shown in Fig. |3(d)| clearly indicates 
the different dynamical time-scales: with a larger infall the 
peaks at core passage are ~ 4.5 Gyr apart, with a smaller 
infall the peaks in mass are separated by ~ 1.7 Gyr. With a 
larger infall there is a marked difference in gas mass before 
and after the peak at core passage suggesting that substan- 
tial stripping occurs during the supersonic phase but that 
this is masked by the density enhancement from the un- 
bound gas. The gas mass variation, with the time axis is 
scaled to the cluster crossing time (Fig. 0, shows that the 
variations in gas mass occur at similar points during the dy- 
namic cycle. The bound gas variations for the above series 
of runs are plotted in Fig. |4(d)| The canonical and smaller 
infall cases have a similar minimum gas content but with 
a smaller infall amplitude the excursions in gas content are 
smaller. 

Analyti c fits t o simulation results have been derived by 
iGaetz et alJ (Il987l) to determine the fraction of replenished 
gas which is retained in a quasi-steady state, but it should 
be noted that the stripping occurring in these simulations 
operat es in a different regime. The results of IGaetz et alJ 
( 1987) are applicable during supersonic regimes only (Mach. 
numbers greater than unity), which applies only briefly in 
these simulations, and in cases where the escape velocity 
from the galaxy (measured at the half mass replenishment 
radius) exceeds the velocity of the galaxy through the ICM. 
In the cases considered here, for the canonical dynamics, 
the escape velocity is comparable wit h the maximum gal axy 
velocity hence the fitting formulae of lGaetz et alJ il987f) do 
not, and would not be expected to, give a good fit to our 
results. 
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Figure 9. Lx:Lb relation for four of the simulated galaxies (less 
massive (triangle), more massive (square), canonical mass with 
lower mass replenishment rate (hexagon) and canonical mass with 
higher mass replenishment rate (star)). Lx is the X-ray luminosity 
in the 0.5—2.0 keV band and Lb is the blue band optical luminosity 
in solar units. The solid line is the observationally determined 
relation of O'Sullivan et al. (2001). Bars in the Lx direction show 
the range in Lx during the simulation (see text for details). 



3.3 Observable halo properties 

The time variation of X-ray luminosity within a radius of 40 
kpc, in the 0.3-8.0 keV band (i.e. a broad band Chandra or 
XMM luminosity) , for runs with varying initial halo extents 
is shown in Fig. |8(a)| These luminosities only include the 
contribution from the gas halo and no point source contri- 
bution. The total X-ray luminosity is seen to broadly follow 
the variations in gas mass within this fixed aperture, i.e. 
peaks at core passage and converges to a common pattern 
of fluctuations when the halo has been stripped. 

The time variation of the galaxy X-ray luminosity for 
different values of the mass replenishment rate is shown in 
Fig. |8(b)| The mass replenishment rate has a large impact 
on the X-ray luminosity; the canonical and higher mass re- 
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Figure 8. Time evolution of X-ray luminosity (Lx) in the 0.3—8.0 keV band for each parameter variation. In each case the canonical 
run is plotted as a solid line, the dashed line represents a reduction in the parameter and the dotted line an increase in the value of the 
parameter. 



plenishment rate runs differ only by a factor of 3 in mass 
replenishment rate but have very different X-ray luminosi- 
ties. 

The time variation of the galaxy X-ray luminosity, for 
different galaxy masses, is shown in Fig. |8(c)] The least mas- 
sive galaxy has a significantly lower X-ray luminosity than 
the canonical case whereas, in contrast, the most massive 
galaxy has a similar luminosity to the canonical galaxy. The 
X-ray luminosity is largely governed by the central regions 
of the gas halo so differences in gas mass may be present 
which do not manifest themselves as significantly different 
X-ray luminosities. 

The time variation of galaxy X-ray luminosity, for dif- 
ferent infall radii, is shown in Fig. |8(d)| The different clus- 
ter crossing times result in different periodicities in the 
Lx variations but it is not possible to distinguish between 
the different dynamical states based on the value of Lx 
alone. An Lx:Lb relation for four of the simulated galaxies 
(lower galaxy mass, higher galaxy mass, canonical mass with 
lower mass replenishment rate, canonical mass with higher 



mass replenishment rate) is presented in Fig. [0] The upper 
and lower bounds indicate the maximum and minimum Lx 
reached during simulation times between 2 and 8 Gyr (to 
avoid any initial transient effects) and the symbols indicate 
the mean Lx measured during the same period. The lower 
mass replenishment rate run is responsible for the lowest 
Lx for a galaxy of canonical mass, and the higher mass re- 
plenishment rate run is responsible for the highest Lx hence 
the inclusion of these two points indicates the range of Lx 
which is induced by parameter variations. The X-ray lumi- 
nosities between 2 and 8 Gyr have been converted from a 
broad band (0.3-8.0 keV) to a soft band (0.5-2.0 keV) to 
enable a direct comparison between X-ray luminosities de- 
termined from ROSAT and the theoretical values derived 
here. A point source contribution, which scales with Lb , has 
been added assuming 



= 29.45 



(17) 



where L ps is the contribution from point sources only 
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JCiotti et alJll99ll : lQ'Sullivan et al.ll200ll) . The X-ray lumi- 
nosities derived for these simulated systems are in agr e ement 
with the observational relation of lO'Sullivan et alJ (2001) 
which is plotted as a solid line in Fig. [§] Variations in Lx 
due to observing a given model galaxy at varying points in 
the dynamical cycle are smaller than those variations in- 
duced by differing model parameters. The largest variation 
for an individual model is for the model with a larger ini- 
tial halo extent which experiences Lx variations of ~ 1.4 
dex, however the effect of parameter variations (particularly 
mass replenishment rate) induces a scatter of ~ 2 dex indi- 
cating that variations in Lx are not solely due to observing 
galaxies at different points in the dynamical cycle but rather 
that there are more fundamental differences causing Lx vari- 
ations. 

Azimuthally averaged X-ray surface brightness profiles 
(0.3-8.0 keV) and temperature profiles, from various times 
during the first and second cluster crossings, are plotted in 
Fig. 1101 The temperature profiles are taken directly from 
the simulation grid so in order for a direct comparison to be 
made between these profiles and observations it would be 
necessary to apply a deprojection analysis to the observed 
profiles. The surface brightness profiles consist, in general, of 
a central peak resulting from a cool, dense core surrounded 
by the cool gas associated with the galactic halo. This is em- 
bedded in the less dense, less bright, hotter ICM which re- 
sults in a str ucture reminiscent of the cold fronts seen in clus - 
ter mergers JVikhlinin et alJl2001al . lMarkevitch et alj|200fj) . 
As the profiles in Fig. 1101 are azimuthally averaged, encom- 
passing regions in front of and behind the galaxy, they lack 
the sharp surface brightness edges seen in Chandra obser- 
vations of merging clusters. When slices along the direction 
of motion, in the upstream direction, are plotted they ex- 
hibit sharp temperature and surface brightness edges (see 
Fig. II IB . The surface brightness slice (Fig. |ll(a)| ) and the 
corresponding temperature slice (Fig. |ll(cJ| l clearly indicate 
the presence of a cold front where a sharp decrease in sur- 
face brightness is coincident with a transition to a higher 
temperature region at log(r/kpc) ~ 1.4. The galaxy ve- 
locity is subsonic at this time so no bow shock is present 
although the gas temperature does rise approaching the dis- 
continuity as would be expected from adiabatic compression 
JVikhlinin et al.ll2001al. §5.3). After 1 .6 Gyr the gala xy mo- 
tion is supersonic and the slices (Fig. |ll(b)| and Fig. |ll(d)] l 
are qualitatively different to the previous subsonic case. 
There is still a cold front present, where a surface bright- 
ness drop coincides with a temperature rise, but at a slightly 
smaller radius (log(r/kpc) ~ 1.2). In addition to the cold 
front there is also a shock front where the temperature and 
surface brightness both rise. The surface brightness enhance- 
ment is significantly greater across the cold front than across 
the shock front which is consistent with Chandra observa- 
tions of Abell 3667 where both a cold front and a shock front 
are detected but the shoc k front results in a sma ller surface 
brightness enhancement JVikhlinin et al.ll2001alh Although 
the size scales and temperatures in this case are different to 
those encountered in cluster mergers the gas dynamics are 
similar to a cluster minor merger. In both the case of a clus- 
ter accreting a group/smaller cluster and the case of a cluster 
accreting a large elliptical galaxy there is a pre-existing halo 
of X-ray emitting gas, confined by a dark matter halo, at a 
temperature below that of the ICM in the accreting system. 




Figure 12. Density (greyscale) and velocity (vectors) around the 
interface between the cool and hot regions. These data are taken 
from the canonical run after 0.8 Gyr. 



IVikhlinin et alJ l|2001rJ) show that conduction, diffu- 
sion and Kelvin-Helmholtz instabilities must be suppressed 
around the Abell 3667 cold front in order for the observed 
sharp discontinuity to be preserved. They propose that a 
magnetic field of order WfiG would be required in a con- 
figuration in which field lines are tangled or run parallel to 
the cold front. Whilst our simulations do not incorporate 
magnetic fields, conduction and diffusion processes will be 
suppressed due to the absence of these terms in the fluid 
equations used. Given the sharp discontinuities observed in 
Chandra data the non-conductive and no n-diffusive fluid 
equatio ns would seem to be appropriate. IVikhlinin et alJ 
( 2001LJ) also suggest that initially tangled magnetic fields 
should be stretched into a morphology where the field lines 
are parallel to the front due to the velocity field of the plasma 
around the interface between the cooler and hotter regions. 
The velocity field seen in our simulations (Fig. I12H clearly 
demonstra tes this and agrees well w ith the flow structure in 
fig. 2(b) of IVikhlinin et alJ feOOlhlh 

Azimuthally averaged surface brightness profiles yield 
useful information regarding the evolution of the gas halo. 
The degree of central concentration of the galaxy halo, out- 
side the central core, may be quantified by fitting a beta 
model to the profile. The surface brightness at a projected 
distance w from the centre of the galaxy is given by 

So 



S(w) = 



\l + (w/r c f 



,3/3- 



(18) 



where r c is a core radius, So is the central surface bright- 
ness and the value of f3 parameterizes the steepness of the 
profile outside the core region. The time variation of the 
beta parameter for fits to the surface brightness distribu- 
tion from all runs is plotted in Fig. 1131 The model fitted 
to the azimuthally averaged surface brightness profile con- 
sists of three components: a beta model to fit the extended 
emission from the galaxy halo, a Gaussian to represent the 
central cool core, and a second order polynomial to account 
for the non-uniform cluster background. The effect of differ- 
ent initial halo extents (Fig. |13(a!J[ l , again shows convergence 
after the initial halo has been stripped. During subsonic mo- 
tion stripping occurs via shear instabilities at the interface 
betweeen the galaxy halo and the ICM, which is less effective 
than stripping in the supersonic phase. This enables replen- 
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(a) Surface brightness profiles during the initial infall, after 
0.8 (solid), 1.6 (dashed) , 2.4 (dotted) Gyr. Surface bright- 
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(b) Surface brightness profiles during the second clus- 
ter crossing, after 3.4 (solid), 4.2 (dashed) , 5.0 (dot- 
ted) Gyr. Surface brightness is measured in units of 
ergs -1 cm -2 arcmin -2 . 
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(c) Temperature profiles during the initial infall, after 0.1 
Gyr (solid), 1.6 Gyr (dashed), 2.4 (dotted) Gyr. 
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(d) Temperature profiles during second crossing, after 3.4 
Gyr (solid), 4.2 Gyr (dashed), 5.0 (dotted) Gyr. 



Figure 10. Azimuthally averaged radial surface brightness profiles and temperature profiles from the canonical run. The centre of the 
galaxy is located at r = 0. Surface brightness S is in units of ergs -1 cm -2 arcmin -2 . 



ished gas to accumulate within the galaxy. Consequently the 
beta value steepens as the regions of the halo where mass 
injection is effective (r < Rh) become brighter. During the 
supersonic phase the beta value drops as gas is stripped 
and the surface brightness profile becomes flatter. This ef- 
fect is not always apparant during the initial crossing but 
can be seen clearly in the third core passage in Fig. |13(a)| 
for example. The effect of mass injection is seen in Fig. |13(b)| 
where higher mass injection rates lead to more rapidly rising 
beta values. This may be attributed to the surface bright- 
ness profile steepening more rapidly due to the increased 
accumulation of gas at radii r < Rh ■ When the infall ampli- 
tude is larger (Fig. |13(d)| , the galaxy spends longer in the 
outer regions of the cluster where the ram pressure is lower, 
so more replenished material can accumulate and the beta 
value grows substantially over an extended period of time. 
The diagnostic power of surface brightness profiles, rel- 



ative to integrated X-ray luminosities, is seen when compar- 
ing the runs with an increased halo mass and an increased 
mass replenishment rate. After 2 Gyr both simulated galax- 
ies have a similar X-ray luminosity but, as Fig. 1141 shows, 
there are marked differences in the surface brightness pro- 
files. The higher mass replenishment rate, which acts only 
within the stellar distribution (log (r/kpc) < 1.3), produces 
a surface brightness profile which is slightly higher in the 
central regions but significantly lower in the outer regions. 
This information regarding the spatial distribution of the 
X-ray emitting gas is lost when considering only integrated 
X-ray luminosities. 
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(a) Surface brightness after 0.8 Gyr in units of 
ergs -1 cm -2 arcmin -2 . The galaxy velocity is subsonic at 
this time. 
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(b) Surface brightness after 1.6 Gyr in units of 
ergs -1 cm -2 arcmin -2 . The galaxy velocity is supersonic at 
this time. 
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(c) Temperature profile after 0.8 Gyr. The galaxy velocity (d) Temperature profile after 1.6 Gyr. The galaxy velocity 

is subsonic at this time. is supersonic at this time. 



Figure 11. On-axis surface brightness and temperature slices from the canonical run. These slices are taken along a line parallel to 
the direction of the galaxy motion through the galaxy centre. The centre of the galaxy is located at r = 0. After 0.8 Gyr the galaxy 
motion is subsonic and a cold front is seen at log (r/kpc) ~ 1.4, after 1.6 Gyr the galaxy motion is supersonic and a bow shock is seen 
at log (r/kpc) ~ 1.6 in addition to a cold front at log (r/kpc) ~ 1.2. 



3.4 Observable wake properties 

In addition to the X-ray emission from the body of the 
galaxy there is also emission from stripped material down- 
stream of the galaxy which may form an observable wake. 
The top row of Fig. ll5l shows X-ray surface brightness maps 
from the canonical simulation during the first, second and 
third crossings of the cluster at the time when the galaxy 
velocity is at a maximum. The highest surface brightness 
feature in all cases is the emission from the body of the 
galaxy and the next brightest is the tail of stripped mate- 
rial. The surface brightness of the tail is much higher during 
the first crossing than in subsequent crossings indicating a 
bias towards observing wakes around galaxies on their initial 
infall into a cluster. The bottom row of Fig. ll5l shows surface 
brightness maps from the run with a larger initial gas halo. 
The wake formed during the first cluster crossing is brighter 



and wider than that from the canonical run suggesting that 
the initial gas content of the galaxy has a major influence 
on whether wakes will be observable. X-ray luminosities, in 
a 0.3-8.0 keV band, for the first cluster crossing are tab- 
ulated in Table [5] These luminosities have been calculated 
as enhancements above the background level which would 
be present at the same distance from the cluster centre but 
away from the galaxy location. The brightness of the wake 
is most strongly affected by the pre-existing halo size; larger 
initial haloes result in more luminous wakes. The mass re- 
plenishment rate does not have a noticable influence on wake 
luminosity during the first infall. Its effect only becomes ap- 
parant at later times when it affects the X-ray luminosity 
from the body of the galaxy. 

The X-ray flux maps presented in Fig. 1151 are valuable 
for determining which features result in the most observable 
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Figure 13. Time variation of beta values from fits to surface brightness profiles. Beta parameterizes the steepness of the surface brightness 
profile outside the central core region. In each case the canonical run is plotted as a solid line, the dashed line represents a reduction in 
the parameter and the dotted line an increase in the value of the parameter. 
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Figure 14. Comparison of the surface brightness profiles for 
runs with an increased halo extent (dotted line) and an in- 
creased mass replenishment rate (solid line) after 2 Gyr in units 
of ergs -1 cm -2 arcmin -2 . 



Table 2. Wake luminosities during the first cluster crossing when 
the galaxy velocity is maximal (i.e. at core passage). 



Run 


log(Lx) 


Canonical 


40.9 


Larger halo 


41.3 


Smaller halo 


40.4 


Larger mass replenishment rate 


40.9 


Smaller mass replenishment rate 


40.9 


Larger galaxy 


40.9 


Smaller galaxy 


40.6 


Larger infall amplitude 


39.9 


Smaller infall amplitude 


40.4 



X-ray emission. However, they should be treated with cau- 
tion as they do not include the effects of detector sensitivity, 
spatial resolution and background levels. Flux maps, such as 
those in Fig. 1151 may be converted to simulated X-ray im- 
ages as shown in Fig. 1161 These simulated images represent 
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(c) Third crossing, t=6.4 Gyr 





(d) First crossing, t=1.3 Gyr 



-200 -150 -100 -50 50 
X (kpc) 



(e) Second crossing, t=3.9 Gyr 



u 

-50 50 100 150 200 
X (kpc) 



(f) Third crossing, t=6.4 Gyr 



Figure 15. Simulated X-ray flux maps from the canonical run (top) and the larger gas halo run (bottom). All images are plotted using 
the same logarithmic greyscale in units of ergs -1 cm -2 arcmin~ 2 . 
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Figure 16. Significance levels in simulated X-ray images derived using the flux maps show in Fig. |15(a)| and Fig. |15(d)| In both cases 
the galaxy is at its first core passage, after 1.3 Gyr of simulated infall. Contours represent significance and are at levels of 1,2,3,4,5 <r. 



50 ksec Chandra observations of the canonical galaxy at first 
core passage (based on the flux map in Fig. |15(a)| l and the 
galaxy with a larger initial halo, also at first core passage 
(based on the flux map in Fig. |15(d)| l. In both cases the 
cluster is placed at a distance of 270 Mpc. An energy band 
of 0.8-2.0 keV is used to increase the contrast between the 
cooler galaxy emission and the hotter cluster background. 
The spatial bins are 5 arcsec square and lkpc = 0.76 arcsec 
at the assumed distance of 270 Mpc. The effects of detector 
background and Poisson noise have been added, a cluster 



background level is then subtracted and the resulting image 
smoothed with a Gaussian filter of 3 pixels FWHM. Even 
at this relatively large distance a tail of material is visible 
in both cases, most prominently in the larger halo case, al- 
though it is affected by noise in the outer regions. These 
simulated images confirm the observability of the structures 
seen in the flux maps: the galaxy halo is the most observ- 
able feature, the stripped wake is the next brightest struc- 
ture (with observability declining with increasing distance 
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from the galaxy) and the bow shock is the least observable 
feature, but may still be detectable. 

A simple statistical measure of the degree of asymme- 
try present in the galaxy's X-ray halo can be formed using 
skewness 

1 \ ' 



7 ; 



Fi,3 - x) 



(19) 
(20) 



where i and j are indices representing cells along the x and y 
axes respectively, Fij is the flux from an individual cell, Ft 
is the total flux from the image, Xij is the location of the cell 
along the a;-axis and x is the x location of the galaxy centre 
along the z-axis. The time variation of the skewness param- 
eter 7, for the larger halo simulation is plotted in Fig. |17(a)1 
This has been calculated using the broad band flux maps. 
Skewness derived from the corresponding simulated images 
is plotted in Fig. |17(b)| for the same run. Fig. |17(a)| rcp- 
resents an idealised case with no noise whereas Fig. |17(b)| 
shows the effects of noise present in real data. When the 
value of skewness is positive the emission is distorted to the 
right and the skewness is negative when the emission is dis- 
torted to the left. Skewness is, in general, a good measure 
of the direction of galaxy motion although when the galaxy 
speed is low, in the outer regions of the cluster, there are 
periods when the skewness does not indicate the correct di- 
rection of motion. Given a sample of galaxies from a cluster 
the skewness parameter will serve as a statistical indicator 
of the direction of motion with respect to the ICM even after 
the first infall phase. A sample drawn from the central re- 
gions of the cluster will preferentially select cases where the 
skewness measure is a reliable direction indicator. In the 
above examples skewness has been calculated in a circular 
region of radius 120 arcsec (~ 90kpc) centred on the galaxy 
to approximately encompass the significant emission from 
any wake formed (c.f. Fig. |15(d)] l. When applying skewness 
to real data the reliability of the measure may be enhanced 
by tuning the size of the region used relative to the extent 
of galactic haloes seen in the observation. 



3.5 Shock properties 

The density and temperature jumps across a shock, such 
as that which precedes the galaxy during supersonic mo- 
tion, can be derived from the Rankine-Hugoniot relations 
l|Landau fc Lifshitzlll95flh . 



p2 

Pi 

T2 

Ti 



2 + ( 7 - 1) Mf 
27M? - 7 + 1 pi 



7 + 1 



f>2 



(21) 



(22) 



where 7 = 5/3 for the monatomic gas considered here. The 
density ratio across the shock front, present after 1.6 Gyr, is 
measured to be 2.0 (see Fig. |18(aJ| which indicates that the 
Mach number of the incident gas is ~ 1.7. The local gas tem- 
perature in the region immediately upstream of the shock is 
~ 2.6 x 10 7 K (see Fig. |l 1 (cTJ| indicating that the gas velocity 
immediately upstream of the shock is ~ 1330 kms -1 , con- 
sistent with the velocity profile in Fig. |18(b)| The predicted 



temperature ratio across the shock is 1.8 consistent with 
that observed in Fig. |ll(d)| The velocity and Mach number 
derived above are local values which pertain to the region 
immediately upstream of the shock. The galaxy is a gravitat- 
ing body and causes gas upstream to accelerate towards the 
shock location hence the velocity of the gas calculated above 
is higher than the global velocity of the galaxy with respect 
to the host cluster. Additionally there is a decrease in the 
local gas temperature during the acceleration towards the 
shock which also causes the Mach number to increase above 
the global value. In the above example the velocity of the 
galaxy with respect to the host cluster is 943 kms -1 but the 
value derived from the density jump is 1330 kms -1 . 



4 STRIPPING IN A LARGER CLUSTER 

The parameter study of a number of scenarios in the same 
modest cluster yields information about how a homogeneous 
sample of wakes may be biased. However, it is also instruc- 
tive to investigate the effect of placing the model galaxy in 
a more substantial cluster where stripping is expected to 
be more effective. Thus we performed a simulation with the 
canonical galaxy placed in a cluster with the cluster tem- 
perature scaled up by a factor of 2 to 5.4 keV. Several pa- 
rameters had to be scaled accordingly, so that the cluster 
remains realistic. The scaled values are shown in Table |3] 
We scale the cluster velocity dispersion a c , the cluster core 
radius r c and the central ICM density nicM as follows. 

ml/! 



to be consistent with eqn. 1141 and 

ml/2 

r c cc T ' 



(23) 



(24) 



assuming that the core radius is a fixed fraction of the virial 
radius and that the virial radius scales with temperature in 
a self similar fashion. The scaling of the central ICM den- 
sity was based on the X-ray emission properties. For a beta 
model surface brightness distribution with f3 — 2/3 the cen- 
tral emissivity So scales as 

So oc j r c (25) 
where the central emissivity jo is assumed to scale as 
jo oc nf CM (26) 
hence 

So cc n\ CM r c . (27) 
From eqn. 1151 

L x oc So 77. (28) 
Combining eg ns, l2"7l and l28l gives 

nf CM oc L x r~ 3 . (29) 

It is assumed that Lx scales as Lx cc T 3 in accordance with 
observationally determined relations llWhite et alJ Il997l) 
rather than Lx oc T 2 which would be expected from self 
similarity. 

The velocity of the galaxy in the hotter cluster is plotted 
as the solid line in Fig. 1191 and for comparision the velocity 
variation for the galaxy in the 2.7 keV cluster is shown by 
the dotted line. The cluster sound speed is shown by the 
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(a) Skewness derived from flux maps 
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Figure 17. Time variation of skewness parameter based on the flux maps from the canonical case. The vertical lines mark regions where 
the galaxy motion is from right to left (R-L) or from left to right (L-R). When the skew is positive the emission is distorted to the right 
and when the skew is negative the distortion is to the left. 
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Figure 18. Density and velocity jumps across the shock after 1.6 Gyr. The both quantities are measured upstream, parallel to the 
direction of galaxy motion along a line through the centre of the galaxy. The galaxy is centred at r = 0. 



Table 3. Parameters of the hotter cluster. 



Parameter 


Value 


Temperature 


5.4 keV 


Velocity dispersion 


928 kms" 1 


Core radius 


566 kpc 


Central ICM density 


1.75 x 10" 3 cm- 3 



solid and dotted horizontal lines for the hotter cluster and 
cooler cluster respectively. The cluster crossing time is the 
same in both cases as galaxy velocity and cluster size both 
scale as T 1/2 . 

The mass of gas gravitationally bound to the galaxy 
(shown in Fig. I20H is significantly less in the hotter cluster 
after the first 2 Gyr. The bound gas content is initially higher 



as the galactic gas halo density must increase to counteract 
the increased ICM pressure confinement. However, this extra 
gas is quickly stripped during the first crossing. The ram 
pressure is higher in a hotter cluster due to the higher ICM 
density and also due to the higher velocities reached. As 
a consequence it would be expected that galaxies residing 
in such an environment would retain less gas than those 
galaxies in cooler clusters where stripping is less effective. 

It is expected that wakes will be more difficult to ob- 
serve in hotter clusters for two reasons. Firstly, the gas con- 
tent of the cluster galaxies should be rapidly depleted by 
stripping, and secondly the background emission from the 
ICM will be brighter than in a cooler system. In order to 
assess the observability of wakes in the hotter cluster we 
have again derived simulated X-ray flux maps and simu- 
lated images. These are presented in Fig. 1211 The flux map 
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X (kpc) 

(a) Simulated X-ray flux map similar to those in Fig. 1151 Note 
the difference in the greyscale levels used here compared to 
those used in Fig, 1151 
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(b) Significance levels in a simulated X-ray image, similar 
those in Fig. 1161 Contours represent significance and are at 
levels of 1,2,3,4,5 a. 



Figure 21. Simulated X-ray flux map and significance map for the hotter cluster run after 1.3 Gyr. 
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Figure 19. Velocity of the galaxy in the hotter cluster (solid 
line) and the 2.7 keV cluster (dotted line). Horizontal lines mark 
the cluster sound speed (solid line is hotter cluster, dotted line is 
2.7 keV cluster). Solid vertical lines mark the transition from the 
subsonic to supersonic regime which is the same in both cases. 



was generated in the same way as those in Fig. 1151 and the 
significance plot was generated in the same way as Fig. 1161 
to enable a direct comparison between the results for the 
hotter and canonical clusters. The gas halo associated with 
the galaxy produces a significant amount of emission, as 
seen in Fig. |21(b)| as the high initial gas content has not 
yet been stripped at this stage. There is a wake of emission 
from stripped material but it is shorter and narrower than 
the previous wakes in Fig. Ilfil indicating that wakes in hotter 
clusters will be less observable than in cooler systems. The 
bow shock may also be detectable. 
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Figure 20. Mass of gas gravitationally bound to the galaxy in 
the canonical (solid line) and hotter cluster (dashed line) runs. 



5 SUMMARY AND DISCUSSION 

We present results from a parameter study involving hydro- 
dynamic simulations of a large elliptical galaxy falling radi- 
ally into an intermediate temperature cluster. The resulting 
stripping causes variations in the gas mass associated with 
the galaxy and in the resulting X-ray luminosity. Variations 
in model parameters (particularly mass replenishment rate) 
are found to induce larger variations in Lx (~ 2 dex) than 
the effects of observing a given model at different points 
during the galaxy's dynamic cycle. 

The luminosity of the X-ray wake during the first clus- 
ter passage is found to be most strongly influenced by the 
extent of the pre-existing gas halo and during subsequent 
passages through the cluster there is not a bright X-ray 
wake. This indicates that samples of wakes in clusters will 
be biased towards large elliptical galaxies, with substantial 
haloes, which are on their first passage through the cluster. 
Wakes surrounding galaxies residing in hotter, more massive 
clusters are, as expected, less easily detected. 
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The infall of a cooler gas halo into a cluster results in 
structures reminiscent of the cold fronts seen in cluster merg- 
ers. The cold front is found to be a more prominent sur- 
face brightness feature than the shock front, present during 
supersonic motion, in agreement with recent Chandra re- 
sults of cluster mergers. The density and temperature ratios 
across the shock front are found to be consistent with values 
from analytic theory but the velocity derived in this way is 
a local value which differs from the global velocity of the 
galaxy with respect to the cluster potential as the galaxy is 
a gravitating body which accelerates the local ICM. Surface 
brightness enhancements and temperature variations aris- 
ing from from the galaxy motion should be detectable with 
present instrumentation such as Chandra and XMM. 
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